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We present a microscopic magnetic model for the spin-liquid candidate volborthite 
Cu3V207(OH)2-2H20. The essentials of this DFT-based model are (i) the orbital ordering of Cu(l) 
3rf3z2_r2 and Cu(2) 3dj.2_y2 (ii) three relevant couplings Jjc, Ji and J2 (iii) the ferromagnetic nature 
of Ji and (iv) frustration governed by the next-nearest-neighbor exchange interaction J2. Our model 
implies magnetism of frustrated coupled chains in contrast to the previously proposed anisotropic 
kagome model. Exact diagonalization studies reveal agreement with experiments. 

PACS numbers: 71.20.Ps, 75.10.Jm, 75.25.Dk, 91.60.Pn 



I. INTRODUCTION 

The search for new magnetic ground states (GS) 
is a major subject in solid state physics. Magnetic 
monopoles in the spin ice system Dy2Ti207 (Ref. Hl-Q), 
the metal-insulator transition in the spin-Peierls com- 
pound TiOCiy and the quantum critical behavior in 
Li2ZrCu04[5|[aare among recent discoveries that demon- 
strate the power of combining precise experimental tech- 
niques with modern theory. However, for a rather large 
number of problems experiment and theory don't keep 
abreast, since it is often tricky to find a real material re- 
alization for a well studied theoretical model. The most 
remarkable example is the concept of a "resonating va- 
lence bond"— — a magnetic GS formed by pairs of cou- 
pled spin-singlets lacking the long range magnetic order 
(LRO). Subsequent studies revealed a fascinating variety 
of disordered GS^^i^ commonly called "spin liquids" in or- 
der to emphasize their dynamic nature, and even raised 
the discussion of their possible applications 

Following the common belief that the spin liquid GS 
may emerge from the interplay of low-dimensionality, 
quantum fluctuations and magnetic frustration, consid- 
erable effort has been spent on the search for spin-1/2 
Heisenberg magnets with kagome geometry. The syn- 
thesis of herbertsmithite Cu3Zn(OH)6Cl2 [HI, the flrst 
inorganic spin-1/2 system with ideal kagome geometry, 
and subsequent studies revealed besides the desired ab- 
sence of magnetic LRO^ (i) intrinsic Cu/Zn structural 
disorder and (ii) the presence of anisotropic interactions 
complicating the spin physicsJ^ The recently synthesized 
kapellasitoi^ was predicted to imply modified kagome 
physics due to an additional relevant couplingJ^ 

Since the search for a system representing the pure 
kagome model is far from being completed, it is natural to 
consider systems with lower symmetry where the distor- 
tion is small enough to keep the essential physics This 
way, the attention has been drawn to the mineral volbor- 
thite Cu3V207(OH)2-2H20, where the Cu sites form a 
slightly distorted kagome networkJi However, the local 
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FIG. 1: (Color online) Top: Cu(l)02 dumbbells (yellow), 
Cu(2)04 plaquettes (dark yellow), V2O7 pyrovanadate groups 
(connected green tetrahedra) and H2O molecules in volbor- 
thite. Only short Cu-O bonds are shown. Bottom: a dis- 
torted kagome layer in the crystal structure of volborthite. 
The magnetically active orbitals and leading exchange cou- 
plings are shown. 



environment of two independent Cu sites is essentially 
different: Cu(l) forms dumbbells of two short Cu-O 
bonds (and four long Cu-O bonds; "2-f4"), while Cu(2) 
resides in a plaquette formed by four short bonds (Fig. [TJ 
top). Recently, DFT studies of CuSb2 06, implying the 
"2-1-4" local environment of Cu atoms, revealed that or- 
bital ordering (00) drastically changes the nature of the 
magnetic coupling from three-dimensional (3D) to one- 
dimensional(lD) The crucial impact of 00 onto mag- 



2 



netism of volborthite will be in the focus of this paper. 

The availability of a pure powder and negligible struc- 
tural disorder in volborthite inspired thorough experi- 
mental studies of its magnetic properties. The magnetic 
GS was recently investigated by ^^V nuclear magnetic 
resonance (NMR),— following earlier NMR,— high- field 
electron spin resonance^i and muon spin relaxations^ 
studies. The GS is characterized by the absence of mag- 
netic LRO, high density of low energy excitations and 
two distinct scales of spin fluctuations. At 4.5 Tcsla, 
volborthite undergoes a transition to another magnetic 
phase. The fingerprint of this transition is a step-like 
feature in the magnetization curve. Similar features are 
observed at 25 and 46 Tesla, hinting at a series of succes- 
sive transitions.— Between 60 and 70 Tesla, the slope 
of magnetization diminishes indicating a possible on- 
set of a magnetization plateau.— Magnetic susceptibility 
measurements yield a broad maximum at temperatures 
much smaller than the Curie- Weiss temperature, indicat- 
ing strong frustration. 

Extensive experimental information stimulated theo- 
retical studies aiming to find a consistent description for 
magnetism of volborthite. Since the pure kagome model 
doesn't account for the experimental data, the studies 
were focused on the GS and thermodynamical proper- 
ties of the anisotropic kagome model (AKM). However, 
attempts to reach consistency by varying the degree of 
anisotropy were not satisfying so far. The most strik- 
ing disagreement is the deviation of the theoretical mag- 
netic susceptibility x even at rather high temperatures 

This disagreement originates from the choice of the 
AKM which was based on geometry only, while the struc- 
tural peculiarities of volborthite were not considered. 
A standard tool to treat such peculiarities properly is 
density functional theory (DFT) calculations that can 
provide a reliable microscopically based modeli^i^ i^^'^^i^^ 
Here, we show that DFT calculations yield an unex- 
pected microscopic magnetic model for volborthite, mov- 
ing away from the kagome model. Moreover, we reveal 
strong similarities to the physics of frustrated coupled 
chains due to 00. Our subsequent simulations of the mi- 
croscopic model evidence an improved agreement with 
the experimental data. 



II. DFT CALCULATIONS 

The DFT calculations have been performed in the local 
density approximation (LDA) using the full potential^^ 
code f plo8 . 65-32<^ For the scalar relativistic calcula- 
tions, the Perdew and Wang parametrization^S of the 
exchange-correlation potential has been used. All cal- 
culations have been performed on well-converged k- 
meshes.— 

The reliability of DFT calculations depends crucially 
on the accuracy of the experimental structural data used 
as input. The chemical composition of volborthite ham- 



j^e A -2 energy (eV) 2 




FIG. 2: (Color online) Density of states (upper panel), the 
band structure (lower panel, center) and the orbital-resolved 
density of states for Cu(l) (left) and Cu(2) (right) of the 
optimized structure {?>d^2_y2 is shown with a dashed line). 

pers structural studies due to the considerable content 
of V and H atoms, which are poor scatterers of neu- 
trons and x-rays, respectively. Therefore, prior to inves- 
tigations of subtle electronic effects, the structural data 
should be addressed. Among several structural data sets 
available we have chosen a structural model (cxp) based 
on joint x-ray and neutron diffraction studies. Although 
such combination improves the reliability of the result- 
ing data, the statistics (number of reflections) is not suf- 
ficient for fully conclusive results. Therefore, we have 
carried out a structural optimization relaxing the atomic 
coordinates and minimizing the forces, since LDA calcu- 
lations for cuprates usually yield accurate and consistent 
resultsJ^ii^— Moreover, to evaluate the influence of the 
different structural models, we perform calculations for 
the experimental as well as for the optimized structure. 

The LDA-optimized crystal structure (opt) yields con- 
siderably lower total energy and atomic forces. Although 
individual bond lengths and angles change up to several 
perccnts (0-H distance increased by ~fO%), the overall 
structural motive and the different local environment of 
Cu(f ) and Cu(2) are inherited from the original model. 

LDA yields a valence band width of 7 eV (Fig. [2]) typ- 
ical for cuprates, and a metallic GS in contrast to the 
green, transparent samples. This well-known problem of 
the LDA originates from the underestimation of strong 
on-site correlations for a Cu ZdP configuration. Never- 
theless, LDA is a reliable tool to evaluate the relevant 
orbitals and couplings For most cuprates. an effective 
one-band model is well justified by a band complex at 
Fermi level (ep) formed by N antibonding bands, where 
A'^ is the number of Cu atoms per cell. In contrast, in 
volborthite {N = 3), six bands in vicinity of ep (Fig. [2]) 
evidence a sizable hybridization of two different 3o? or- 
bitals at each Cu site that need to be included into the 
modeling. 
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The relevant Cu 3d orbitals are revealed by project- 
ing the density of states (DOS) onto local orbitals. The 
resulting orbital-resolved DOS is shown in Fig. [2] For 
both Cu(l) and Cu(2), the 8^3.2 _,y2 and 3c?322_^2 states 
are relevant and hybridized with each other. To evaluate 
the relevant couplings, we consider two orbitals (Sd^^^-j/^ 
and 3c?3j,2_r2) per Cu atom and fit the six bands using 
the Wannier functions (WFs) tcchnique<^i^ 

Prior to evaluation of the relevant couplings, the cor- 
rect orbital GS should be found. LDA yields an es- 
sentially different filling of the orbitals: for Cu(l) the 
Sd^z^-r^ is close to half-filling and the 3c?a;2_j^2 is almost 
filled (Fig. [H left), while for Cu(2) it is the other way 
round (Fig. [2l right). However, the closer proximity to 
half-filling in the LDA picture does not necessarily pro- 
vide the correct answer, as revealed for the related sys- 
tem CuSb206'^ Thus, we cross-check the LDA result by 
LSDA+U calculations. In agreement with the LDA, the 
latter yield the magnetically active Cu(l) 3d^z^_j.2 and 
Cu(2) 3d^2_y2 (sec details below). 

The relevant transfer integrals can be extracted from 
the WFs considering the hoppings between the GS or- 
bitals (Cu(l) 3(i322_r2 and Cu(2) 3dx2-y2). The leading 
terms are ti and t\c, which coincide with two NN cou- 
plings in the AKM (Fig. [U bottom). Surprisingly, the 
next-nearest-neighbor (NNN) coupling t2 is also sizable, 
while other couplings are considerably smaller. A small 
value of the inter-layer coupling supports the 2D nature 
of magnetism. 

The correct description of the orbital GS requires an 
appropriate description of correlations in the Cu 3d shell, 
which can be treated in a mean-field way using the 
LSDA-|-[/ scheme. By stabilizing solutions comprising 
different orbital occupations and a subsequent compar- 
ison of their total energies, we evaluate the orbital GS. 
The separation between the orbital GS and the lowest 
lying excited orbital state {3d^2_y2 for Cu(l) and Cu(2)) 
exceeds 500 meV (■^GOOO K), almost two orders of magni- 
tude larger than the magnetic exchange (-^100 K). There- 
fore, orbital and spin degrees of freedom are mostly de- 
coupled and can be analyzed separately. 

The leading exchange integrals Ji, J2 and Jjc are ob- 
tained mapping the results of LSDA+J7 total energy cal- 
culations onto a Heisenberg modeli^ A careful analysis of 
the results shows that the individual values of exchange 
integrals are sensitive to (i) the structural model, (ii) the 
Coulomb repulsion Ud and (iii) the double-counting cor- 
rection (DCC) schemei^ The crucial infiuence of these 
parameters is visualized in Fig. [3l where the results for 
the experimental and optimized structures are shown. 
For each structural model, we use the limiting cases for 
the DCC — around-the-mean-field (AMF) and the fully 
localized limit (FLL)2& and vary Ud within a reason- 
able range^i Depending on these parameters, we obtain 
Ji = -80±10 K, J2 = 35±15 K and Jic = I00±60 K for 
the experimental and Ji = — 65±15 K, J2 = 45±I5 K 
and Jic = 100±60 K for the optimized lattice. 

The substantially ferromagnetic nature of Ji, in ac- 
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FIG. 3: (Color online) Left: exchange integrals as a func- 
tion of the structural model (opt or exp), the LSDA-f f/ DCC 
scheme (AMF or FLL) and Ud (min Ud, max Ud) on the phase 
diagram of the Ji- J2- Jic model. The white and gray fields cor- 
respond to the singlet and the ferrimagnetic phases, respec- 
tively. The shaded areas depict possible values of exchange 
couplings for the both structural models. Our ED results yield 
the best agreement for Jic/|Ji| = 2, J2/IJ1I = 1.1, depicted 
by a cross. Right: correlation functions along (squares) and 
between (circles) the J1-J2 chains deviate significantly from 
the kagome model (bold gray Une). 

cord with Goodenough-Kanamori- Anderson rules, and 
the relatively small uncertainties of its strength disre- 
garding the parameters used give strong evidence that 
the pure kagome model is inappropriate for volborthite. 
The AKM can be ruled out since Ji and Jic support 
each other and do not give rise to frustration. Our mi- 
croscopic insight evidences that an essentially different 
model with frustration governed by NNN exchange J2 
(competing with both Ji and Jc) should be used for 
volborthite. Despite sizable scattering of the J values, 
important general trends can be established. First, the 
optimized structure has an enhanced J2 / 1 Ji | ratio com- 
pared to the experimental structure. Second, FLL yields 
considerably smaller Jc and somewhat larger values for 
J2 than AMF. 

Based on DFT calculations, we obtain a microscopic 
magnetic model and determine the parameters Ji , J2 and 
Jic. Although wc find the relevant region in the phase 
space, the complexity of volborthite impedes a more ac- 
curate estimate of individual exchange integrals, espe- 
cially Jic. In this case, refining the parameters by nu- 
merical simulations of measured physical properties and 
subsequent comparison to experimental data is an appro- 
priate way towards a deeper understanding. 



III. EXACT DIAGONALIZATION 

To realize a guided search for a consistent set of ex- 
change integrals, the GS of the Ji-J2-Jic model should 
be investigated. We explore the phase space by consid- 
ering a Heisenberg model with Ji < 0, J2 > 0, Jc > 0. 
On a classical level, we find two GS: a ferrimagnetic (fM) 
phase with magnetization m = 1/3 and an incommen- 
surate m = helical (H) phase with spiral correlations 
along J1-J2 frustrated chains, similar to those of edge 



4 




FIG. 4: (Color online) Top: fits to the experimental x{T)M 
The solution of the Ji-J2-J\c model yields an improved de- 
scription down to 50 K compared to the kagome model (bold 
gray line). Bottom: magnetization curves (A'^=36 sites) for 
different solutions of the Ji-J2-Jic model in comparison to 
the kagome model. 

shared quasi ID cuprates (see e.g. Ref. [^. The tran- 
sition from the fM-phase to the H-phase is driven by 
the frustrating NNN in-chain coupling J2 and occurs at 
jciass ^ I _^ j^jg^ discuss the GS of the quan- 
tum model we use Lanczos exact diagonalization of finite 
lattices up to iV = 36 sitesi^i^ For the quantum model, 
the fM state competes with a singlet GS with m = 0, and 
the transition is given by Jl"""* = 0.304|Ji| +0.200 J^^^a 
The transition line together with the DFT-derived ex- 
change integrals are plotted in Fig. [3] (left). Since the 
experiments evidence zero magnetization of the GS, the 
fM solutions can be ruled out and the analysis can be 
restricted to the singlet GS. 

To understand the nature of the GS, we consider spin 
correlations as a sensitive probe for magnetic ordering. 
While the correlations between the chains are similar to 
the standard kagome Heisenberg antifcrromagnet, they 
are completely different along the J1-J2 frustrated chains 
(Fig. m right). These in-chain correlations fit to a spiral 
state with a pitch angle very close to the classical model. 
Hence, our data suggest well pronounced in-chain spi- 
ral correlations together with weaker inter-chain correla- 
tions. We mention, however, that these statements are 
restricted to short-range correlations. 

Since the magnetic correlations along the chains are 
strongest, one could argue that the model exhibits an 
effectively ID low-temperature physics as has been dis- 
cussed previously for other 2D models such as the 
crosscd-chain model^ anisotropic triangular lattice;^ as 
well as for modified kagome compoundsJ^ However, this 



issue as well as a conclusive answer to the question of 
helical LRO need further investigation. 

To comprise the experimental magnetization curve, we 
add the magnetic field term to the Heisenberg Hamilto- 
nian and simulate the m{h) dependence. For the bound- 
ary of the fM and singlet GS, wc find a wide 1/3 mag- 
netization plateau starting at hd = 0. However, the 
modification of the exchange parameters, in particular 
increasing of J2 and decreasing of | Ji | and J\c , according 
to the limits set by DFT calculations leads to a signifi- 
cant increase of hd and to a drastic diminishing of the 
plateau width. Close to the DFT-boundary (Jic/| Ji| = 2, 
J2/I Jil 1.1, = 100 K), we obtain hd = 22 T, which 
is still smaller than the experimentally observed value. 
We should note that this deviation originates from the 
minimalistic character of the model and considerable fi- 
nite size effects. Nevertheless, a slightly modified ra- 
tio J2/IJ1I = 1.6 yields hd = 55 T (Fig. HI bottom) 
in excellent agreement with the experiment. It should 
be mentioned that the nature of spin correlations in the 
1 /3-plateau phase is substantially different compared to 
the kagome model.— Small magnetization jumps seen 
experimentally^ can not be resolved with present lattice 
sizes and might be related to anisotropic exchange. 

We also calculate the temperature dependence of mag- 
netic susceptibility x(T) using two different lattices^ up 
to TV = 24 sites. We obtain a good fit down to 50 K 
(Fig.H top)^ the resulting g = 2.16 and Jic = 100.5 K 
are in excellent agreement with experiments^! and our 
estimates from DFT. 



IV. SUMMARY 

To summarize, we suggest a new magnetic model for 
volborthite: assuming that DFT calculations arc appli- 
cable to volborthite as they arc for a plethora of com- 
pounds, the kagome model can be safely ruled out. In- 
stead, the magnetism of volborthite is accounted for 
by a J1-J2-J1C model reminiscent of coupled frustrated 
chains. For the proposed model, the orbital order of 
Cu(l) 3d:iz'2-r^ and Cu(2) 3dx2-y2 orbitals is crucial. 

We suggest new experiments to challenge our model: 
resonant x-ray scattering measurements measurements to 
study orbital effects and measurements in high magnetic 
fields (>70 Tcsla) to get an access to the magnetization 
plateau. Additional investigations of the Ji-J2-J\c model 
itself by alternative simulation methods are desirable to 
clarify the influence of finite size effects, intrinsic for exact 
diagonalization. 
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FIG. SI: Finite lattices used for exact diagonalization on A'^ = 18 and N = 24: sites. 



TABLE SI: Basis set used in fplo8.65-32. 



atom states 

core valence second valence polarization 

Cu Is 2s 2p 3s 3p 4s 3d 5s 4d 4p 

V Is 2s 2p 3s 3p 4s 3d 5s 4d 4p 

O Is 2s 2p 3s 3p 3d 

H Is 2s 2p 
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FIG. S2: Comparison of the band structures calculated within the local density approximation (LDA) and the general gradient 
approximation (GGA) for the experimental ("exp", top panel) and optimized ("opt", bottom panel) structures. 



TABLE S2: Cells used for LDA and GGA (DFT) as well as spin-polarized LSDA+C/ and GGA+U (DFT+C/) calculations. For 
each cell, the total number (tot.) as well as the number of symmetrically inequivalent (ineq.) Cu atoms are provided. Basis 
vectors are given in terms of the unit cell vectors a, b and c. PUC — primitive unit cell, SC — supercell. 



cell 


Sp. gr. 


functional 


Cu 


atoms 




basis vectors 




fc-mesh 








tot. 


ineq. 


0.5(a + 6) 


0.5(a - b) 






PUC 


P2/m 


DFT 


3 


2 


c 


12x12x12 


SCI 


PI 


DFT 
DFT+f/ 


3 


3 


0.5(a + 6) 


0.5(a - b) 


c 


8x8x8 
4x4x4 


SC2 


Pi 


DFT 
DFT+f/ 


6 


6 


0.5(3 + b) 


(a-b) 


c 


6x3x6 
3x2x3 


SC3 


PI 


DFT 

DFT+C/ 


6 


6 


a 


b 


c 


4x7x5 
2x4x3 
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FIG. S3: The Wannier functions fit (WF) and the tight-binding fit (TB) together with the LDA band structure for the 
experimental ("exp", top panel) and optimized ("opt", bottom panel) structures. 



TABLE S3: Transfer integrals t as a function of the structural model (exp or opt) and the exchange-correlation potential 
(LDA or GGA), calculated using the Wannier functions technique. The numbers given correspond to hoppings between the 
magnetically active orbitals (Cul 3(1^^2 and Cu2 3d^2_y2). 



exp opt 

LDA GGA LDA GGA 

ti, meV 91 93 117 119 

t2, meV 59 57 64 62 

tic, meV 156 156 155 157 



